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L  INTRODUCTION 


This  is  a  project  to  provide  mission  support  for  the  Communication/Navigation  Outage  Forecast 
System  (C/NOFS)  under  BAA  VS-03-G1  during  its  first  four  years  of  operation.  Cornell  will 
support  the  mission  with  ground-based  radar  observations  of  background  ionospheric  parameters 
and  of  equatorial  spread  F  (ESF)  events  from  the  Jicamarca  Radio  Observatory  near*  Lima,  Peru. 
Jicamarca  is  capable  of  measuring  ionospheric  electric  field  and  conductivity  profiles  from  the 
valley  region  well  into  the  topside  ionosphere.  Jicamarca  can  also  make  very  detailed  observations 
of  ionospheric  irregularities  associated  with  ESF.  Three  main  tasks  will  be  supported  under  this 
project:  (1)  in-flight  calibration  of  the  electric  field /drift  meter  sensors  on  board  the  spacecraft, 
(2)  provision  of  Jicamarca  radar  support  for  experimental  campaigns,  both  before  the  launch  and 
during  the  satellite  operations,  and  (3)  analysis  and  interpretation  of  the  data  obtained  during 
these  campaigns  in  fight  of  the  C/NOFS  mission  goals. 

2.  METHODS  AND  PROCEDURES 

The  C/NOFS  data  workshop,  co-chaircd  by  Hyscll,  took  place  in  January,  2005,  in  Estes  Park, 
Colorado.  At  the  meeting,  a  new  theory  regarding  ionospheric  preconditioning  and  the  role  of  shear* 
flow  in  destabilizing  the  postsunset  equatorial  ionosphere  was  presented.  The  theory,  established  in 
part  using  data  from  C/NOFS  prc-laxmch  preparations  at  Jicamarca,  has  important  implications 
for  forecasting  ESF  and  for  the  conduct  of  the  C/NOFS  mission  and  is  described  in  two  publications 
supported  by  this  award: 

Hyscll,  D.  L.,  and  E.  Kndeki,  Collisional  Shear  instability  in  the  equatorial  F  region  iono¬ 
sphere,  J.  Geophys .  Res 109 ,  10J029//2004JA0106362004,  2004. 

Hyscll,  D.  L.,  E.  Kndeki,  and  J.  L.  Chau,  Possible  preconditioning  of  the  equatorial  ionosphere 

by  shear  flow  leading  to  spread  F,  Annales  Geophysicae,  in  press,  2005. 

The  theory  is  summarized  below. 

Fig.  1  shows  examples  of  ionospheric  structure  and  circulation  from  paired  events  observed  at 
Jicamarca  in  November  or  December  2001,  2002,  and  2003.  The  top  panel  shown  for  each  event 
represents  relative  backseat  ter  power  in  grayscale  logarithmic  format.  Regions  of  intense  backscat- 
tcr  signify  coherent  scatter  from  plasma  irregularities  and  arc  plotted  using  different  amplitude 
scaling.  The  middle  panel  of  each  event  in  Fig.  1  shows  vertical  velocity  derived  from  the  average 
Doppler  shift  of  Jicain area's  dual  beams,  which  arc  directed  approximately  3  degrees  cast  and  west 
of  magnetic  zenith  for  this  experiment.  The  bottom  panel  shows  zonal  velocity  derived  from  the 
difference  of  the  dual  beam  Doppler  shifts.  The  three  examples  on  the  left  side  of  Fig.  1  repre¬ 
sent  events  when  ionospheric  irregularities  were  confined  to  narrow  layers  in  the  bottomsidc.  In 
contrast,  the  three  examples  on  the  right  portray  occurrences  of  ‘Tull-blown”  equatorial  spread 
F.  Forecasting  topside  spread  F  events  like  these  on  the  basis  of  available  data  is  a  goal  of  the 
C/NOFS  mission.  Forecast  strategies  have  often  involved  applying  a  threshold  condition  to  the 
intensity  of  the  prcrcvcrsal  enhancement  of  the  vertical  drift.  Cursory  examination  of  the  data  here 
alone  suggest  that  such  a  strategy  should  have  merit  but  would  be  fallible  and  could  not  be  applied 
before  about  2330  UT  (1830  LT)  or  only  about  an  hour  before  the  onset  of  instability. 
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Figure  1:  Six  examples  of  ionospheric  conditions  around  twilight.  Each  example  shows  the  backseat- 
ter  power  (top  panel),  vertical  plasma  drift  (middle  panel),  and  zonal  plasma  drift  (lower  panel), 
plotted  against  the  scales  shown.  Symbols  plotted  in  the  middle  panel  represent  the  vertical  drift 
velocity  at  450  km  altitude  for  reference.  The  left  and  right  columns  portray  events  when  topside 
spread  F  did  not  and  did  occur.  respectively.  Note  that  UT  =  LT  +  5h. 
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Figure  2:  Range  time  intensity  (RTI)  image  of  a  topside  spread  F  event.  Only  a  bottom-type  layer 
was  present  prior  to  about  2100  LT. 


Shear  flow  is  evident  in  all  of  the  events  in  Fig.  1  and  is  most  obvious  after  about  0030  UT 
(1930  LT)  when  backseat  ter  from  bottom- type  layers  points  to  rapid  westward  motion  immediately 
below  regions  of  rapid  eastward  motion.  It  might  seem  fortuitous  that  the  hot  tom- type  layers  exist 
to  provide  a  radar  target,  since  the  valley  region  plasma  is  insufficiently  dense  to  be  monitored 
using  incoherent  scatter.  However,  [4]  pointed  out  that,  assuming  the  thermospheric  winds  arc 
eastward  throughout  the  F  region,  the  bottom- type  layers  occupy  strata  where  the  plasma  drifts 
arc  strongly  retrograde.  Given  that  significant  horizontal  density  gradients  arc  also  present  in 
these  strata  around  sunset  due  to  photochemical  and  dynamical  effects,  the  strata  would  be  prone 
to  horizontal  wind-driven  interchange  instabilities.  They  hypothesized  that  such  instabilities  arc 
responsible  for  the  bottom-type  scattering  layers. 

Evidence  that  wind-driven  instability  is  at  work  in  the  bottomsidc  was  provided  by  [1],  who  used 
aperture  synthesis  imaging  techniques  to  analyze  the  structure  of  the  primary  waves  in  bottom- 
type  scattering  layers.  The  primary  waves  they  observed  appeared  to  grow  by  advcction  rather 
than  convection,  consistent  with  the  wind  driven  interchange  instability  theory.  Moreover,  [1] 
foimd  that  the  primary  plasma  waves  in  bottom- type  layers  sometimes  cluster  together  in  regular, 
wavelike  patches.  Clustering  took  place  only  in  those  events  when  topside  spread  F  occurred 
later.  The}'  associated  the  clustering  with  a  horizontal  wave  in  the  background  plasma  density 
that  would  present  alternately  stable  and  unstable  horizontal  gradients  for  wind-driven  interchange 
instability.  Fig.  3  shows  a  radar'  inrage  derived  from  the  bottom- type  layer  in  Fig.  2  illustrating 
such  clustering.  Finally,  [1]  concluded  that  the  clustered  or  patchy  bottom-type  layers  were  telltale 
of  dccakilomctric  waves  in  the  bottomsidc  F  region  that  could  serve  to  precondition  or  seed  the 
ionosphere  for  gravitational  interchange  instabilities  (Rayleigh  Taylor)  to  follow.  This  hypothesis 
was  supported  by  the  fact  that  the  intermediate-scale  plasma  irregularities  that  formed  later  in  the 
spread  F  events  shared  the  same  dccakilomctric  wavelength.  This  suggests  a  forecast  strategy  for 
spread  F  based  on  detecting  clustering  in  bottom- type  layers  innncdiatcly  after  sunset. 

Note  that  the  30  km  wavelength  of  the  clustering  is  roughly  comparable  to  the  vertical  scale 
length  of  the  plasma  shear  flow.  Recently,  [2]  examined  the  possibility  that  shear  flow  itself  may 
cause  the  bottomsidc  F  region  to  become  mis  table  and  produce  the  precursor  waves  inferred  from 
radar  imaging.  If  so,  then  shear  flow  could  precondition  the  ionosphere  for  full-blown  spread  F. 
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Figure  3:  Radar  image  of  plasma  irregularities  within  the  bottom-type  layer  seen  in  the  previous 
figure.  The  irregularities  arc  shown  at  the  moment  they  first  drifted  westward  into  the  radar 
illuminated  volume. 

[2]  addressed  the  stability  of  the  F  region  ionosphere  in  a  sheared  flow  with  an  eigenvalue 
analysis.  Following  [3],  [2]  neglected  Hall  and  parallel  currents  but  allowed  for  polarization  currents 
in  their  two-dimensional  analysis: 


a) 


j 


The  quasincutrality  condition,  together  with  the  ion  continuity  equation,  formed  the  basis  of  their 
model.  Linearization  of  the  model  was  performed  according  to 


n(xyz)  =  n0(z)  +  n\{z)e%(k*x 
<t>{x,z)  =  Mz)  + 


(2) 

(3) 


where  terms  with  subscripts  0  and  1  represent  zero-  and  first-order  quantities,  respectively,  and 
where  the  Cartesian  coordinates  x ,  y.  and  z  represented  the  eastward,  northward,  and  upward 
directions  at  the  magnetic  equator.  The  linearized  model  equations  took  the  form: 


(4) 


+tT  {Bu m(u  -  va)ni) 
az 


(5) 
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where  uj  is  the  complex  frequency  and  serves  as  the  eigenvalue,  kx  is  the  assumed  horizontal 
wavenumber,  v0  is  the  background  horizontal  plasma  flow  speed  arising  from  gradients  in  the 
zero-order  electrostatic  potential,  and  vm  is  the  ion-neutral  collision  frequency. 

This  model  was  solved  computationally  for  equilibrium  velocity,  collision  frequency,  and  density 
profiles  specified  by: 


v0(z) 
Vm  {z) 

nQ(z) 


V0Tznh(z/L) 
v0  —  zqTanhl^/L) 
Vo/V]n{z) 


N0u 


U  —  Vo  (z) 


and  shown  to  yield  rapidly  growing  solutions  even  in  the  collisional  regime  (^m  >*  |o;|). 

An  approximate,  closed- form  expression  for  the  growth  rate  highlighting  its  dependence  on 
background  ionospheric  parameters  is  derived  below.  In  our  analysis,  we  focus  on  the  collisional 
branch  of  the  instability  thought  to  operate  in  the  bottomsidc  equatorial  F  region  by  applying  the 
limit  i  oo  in  (4)  and  (5),  leading  to: 


_d_ 

dz 


(  d<f>  A 

{^n°  dFj  ~ 


-kx 


dz 


V\r 


V\j\tl0kx(j)Y 

u  —  v0  dnQ 
1uj  —  kxv o  dz 


fa 


We  next  multiply  the  equation  by  <f>\  and  integrate  over  altitudes  where  the  perturbed  potential 
exists: 


-  J  dz  iflDn0  (|$|2  + 

=  kx  f  dz  wm  -  n0'<f>i4>\'  (6) 

J  u>-  kxv0 


where  the  prime  notation  refers  to  differentiation  with  respect  to  z  and  where  integration  by  parts 
has  been  employed.  Since  the  left  side  of  (6)  is  real,  the  frequency  uj  =  ur  -j-  17  can  only  be  complex 
valued  if  the  eigenfunction  <f>[  is  also  somewhere  complex. 

The  numerical  calculations  of  [2]  suggested  that  in  the  vicinity  of  the  shear  node,  where  v” 
vanishes  and  about  which  <t>  is  localized,  the  complex  potential  can  be  approximated  by  (j)(z )  ^ 
<f)0  cxp(ikzz).  Utilizing  this  approximation  and  neglecting  ujr  —  kxvQ  by  comparison  to  7  in  this 
region,  an  expression  for  the  growth  rate  7  can  be  derived  from  the  real  part  of  (C): 

kxkz  f  dz  (ti  —  Vo)n0f\<f>y\2 
7  ~  k2  f  d.z  v-mn0\<j>i\2 

kxkjjyjnju  -  t^nc/)  .  . 

*>inn0)  1  1 

where  the  averaging  is  weighted  by  the  norm  of  the  potential,  a  function  that  generally  peaks  about 
one  scale  height  below  the  shear  node  and  extends  for  several  scale  heights  above  and  below  it. 

Given  mode  shapes  and  values  for  kz  derived  from  a  complete  boundary  value  analysis,  (7) 
gives  growth  rate  estimates  in  good  agreement  with  the  eigenvalues  of  the  full  model  expressed 
by  (4)  and  (5).  The  final  expression  shows  that  shear  instability  is  likely  to  occur  in  layers  with 


5 


rapid  retrograde  plasma  motion  collocated  with  a  steep  vertical  plasma  density  gradient.  This  is  a 
testable  hypothesis. 

The  task  remaining  is  to  test  and  validate  the  shear  instability  model  using  existing  ground  based 
data  and,  when  available,  in  situ  observations  from  C/NOFS.  We  can  then  assess  the  importance 
of  shear  flow  in  determining  the  stability  of  the  postsunset  ionosphere.  The  capabilities  of  the 
Jicamarca  radar,  along  with  samples  of  our  data  products  and  their  relevance  to  this  model,  were 
presented  to  the  spring  C/NOFS  science  team  at  Hanscom  and  the  June  CEDAR  meeting  in  Santa  * 

Fc:  where  Jicamarca  scheduling,  calibration  and  validation  issues,  and  larger  issues  pertaining:  to 
ionospheric  stability  and  scintillation  forecasts  were  also  discussed.  We  arc  presently  upgrading 
the  radar  mode  with  winch  we  will  support  the  satellite  after  launch  and  planning  our  modeling 
strategy. 
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